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Abstract 

Many cell types exhibit oscillatory activity, such as repetitive action potential fir- 
ing due to the Hodgkin-Huxley dynamics of ion channels in the cell membrane or 
reveal intracellular inositol triphosphate (IP3) mediated calcium oscillations (CaOs) 
by calcium-induced calcium release channels (IPs-rcccptor) in the membrane of the 
endoplasmic reticulum (ER) . The dynamics of the excitable membrane and that of the 
IPs-mediated CaOs have been the subject of many studies. However, the interaction 
between the excitable cell membrane and IPs-mediated CaOs, which are coupled by 
cytosolic calcium which affects the dynamics of both, has not been studied. This study 
for the first time applied stability analysis to investigate the dynamic behavior of a 
model, which includes both an excitable membrane and an intracellular IPs-mediated 
calcium oscillator. Taking the IPs concentration as a control parameter, the model 
exhibits a novel rich spectrum of stable and unstable states with hysteresis. The four 
stable states of the model correspond in detail to previously reported growth-state 
dependent states of the membrane potential of normal rat kidney fibroblasts in cell 
culture. The hysteresis is most pronounced for experimentally observed parameter val- 
ues of the model, suggesting a functional importance of hysteresis. This study shows 
that the four growth-dependent cell states may not reflect the behavior of cells that 
have differentiated into different cell types with different properties, but simply reflect 
four different states of a single cell type, that is characterized by a single model. 
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Complexity and multiple transitions among beliaviorial states are ubiquitous in biologi- 
cal systems (1, 3)- In physics instabilities and hysteresis are well known to play an important 



role in collective properties and have been studied since many years 



). Recently, 



multi-stability with hysteresis has also awakened a large interest in biological systems 
Instabilities, for instance, are crucial for efficient information processing in the brain, such as 
in odor encoding Moreover, unstable dynamic attractors have been demonstrated in 



cortical networks, with critical relevance to working memory and attention ()10llllL 11211 . In a 
wide sense, multistable systems allow changes among different stable solutions where the sys- 
tem takes advantage of instabilities as gateways to switch between different stable branches 
(Q). Bistability driven by instabilities prevents the system from reaching intermediate states, 
e.g. partial mitosis. Hysteresis prevents the system from changing its state when parameter 
values, that characterize the system, vary. This is of relevance, for instance, in cell mitosis. 
Once initiated, mitosis should not be terminated before completion (jl^i . Thus, hysteresis 
may lock the cell into a fixed state, preventing it from sliding back to another state (jl^ . 

At the level of cell networks, multistability, and in particular bistability, plays an im- 
portant role in cell signaling as well (Q, Q). For example, communication between neurons 
takes place at synaptic contacts, where arrival of an action potential stimulates release of 
a neurotransmitter, thus affecting the post-synaptic potential of the target cell. Typically, 
each cell receives input from thousands of other cells mediated by different neurotransmit- 
ters, which modify the post-synaptic potential by excitation or inhibition at different time 
scales jiyi . This information at the cell membrane may be transferred to the cell nucleus 
by so-called second messengers to affect the nucleus in controlling DNA-expression, protein 
synthesis, mitosis, etc. Calcium is one such second messenger and calcium oscillations have 
been reported over a wide range of frequencies with a chaotic or regular pattern ( isi. 

In many biological systems, cells display spontaneous calcium oscillations (CaOs) and 
repetitive action-potential firing. These phenomena have been described separately by mod- 
els for intracellular inositol trisphosphate (IP3)-mediated CaOs jiol l2fli and for plasma 
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membrane excitability (Q). We have recently presented a single-cell model that combines 
an excitable membrane with an IPs-mediated intracellular calcium oscillator The IPs- 
receptor is described as an endoplasmic reticulum (ER) calcium channel with open and close 
probabilities that depend on the cytoplasmic concentrations of calcium ([Ca^^]) and IP3 
([IP3]). An essential component of this model relates to store-operated calcium channels in 
the plasma membrane. Since it is not known whether multiple types of store-operated cal- 
cium channels are involved in normal rat kidney (NRK) fibroblasts, we will use the general 
terminology of store-dependent calcium (SDC) channels. 

NRK fibroblasts in cell culture exhibit growth-state dependent changes in their electro- 
physiological behavior (23). Subconfiuent-grown serum-deprived quiescent cells exhibit a 
stable resting membrane potential near -70 mV ("resting state"). Upon subsequent treat- 
ment with epidermal growth factor the cells re-enter the cell cycle, undergo density-dependent 
growth-arrest (contact inhibition) at confluency and spontaneously fire action potentials as- 
sociated with intracellular calcium oscillations ("AP-firing state"). Subsequent addition of 
retinoic acid or transforming growth factor (TGF)/3 to the contact inhibited cells causes the 
cells to become phenotypically transformed and to depolarize the cell to approximately -20 
mV ("depolarized state"). This depolarization has been shown to be caused by an eleva- 
tion of the concentration of prostaglandin (PG)F2q secreted by the unrestricted proliferating 
transformed cells. Washout of the medium conditioned by the transformed cells by perfusion 
with fresh serum-free medium causes the cells to slowly repolarize, and, preceded by a short 
period of fast small-amplitude spiking of their membrane potential ("fast oscillating state"), 
to regain spontaneous repetitive action potential firing activity ("AP-firing state") similar 
to that of the contact inhibited cells. These phenomena have been described in detail (Q) 
and are very similar to the behavior of other cell types with calcium oscillations and action 

n n 

potential firing, such as interstitial cells of Cajal (|2J) and hepatocytes ()25[). 

In this study we have analyzed the model reported in ()22l ). This model, which is shown 



schematically in Fig. ^ illustrates the basic characteristics of NRK fibroblasts. It repro- 
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duces, on the basis of single-cell data (j2 



^ 2^ , the dynamics of both the plasma membrane 
excitability and that of the intracellular calcium oscillator. We have recently shown that 
(PG)F2a dose-dependently induces IPs-dependent intracellular calcium oscillations in NRK 



fibroblasts (j27|). Since the growth-state dependent modulation of the membrane potential 
of NRK fibroblasts is related to the concentration of (PG)F2a in their culture medium (Q) 
and since this prostaglandin dose-dependently increases [IP3], we took [IP3] as a control 
parameter to analyze the stability of the single-cell model. 

The stability analysis shows how coupling of an excitable membrane with an intracellular 
calcium oscillator leads to a rich behavior of a cell with multiple stable and unstable states 
with hysteresis. We show that the growth-state dependent modulations of the membrane 
potential of NRK fibroblasts in cell culture described above can be understood as the stable 
states of the single-cell model with membrane excitability and calcium oscillations of these 
cells. The stable states of the model reproduce the four growth-dependent states of NRK 
cells, corresponding to the resting state at —70 mV, the AP-firing state for spontaneous 
action potential firing, the depolarized state at —20 mV and the fast oscillating state with 
small-amplitude spiking around —20 mV. Therefore, the four growth-dependent states of 
NRK fibroblasts may not refiect the behavior of cells that have differentiated into different 
cell types with different properties, but reflect four different states of a single cell type, that 
is characterized by a single model. 

Model description 

The dynamics of NRK cell membrane excitability is given by a set of equations which describe 
the active and passive ion transport systems in the plasma membrane and the endoplasmic 



reticulum, as illustrated in Fig. [T] (see (j21 



22) for a detailed description). The change in the 



membrane potential as a function of time due to the currents through inwardly rectifying 
potassium channels (Ixir), L-type Ca-channels (Icai), Ca-dependent Cl-channels {Ici{Ca)), 
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leak channels (Iik), and SDC-cliannels (Isdc) is given by 



dVr 



Cm-^^ — —{I Kir + hk + ^CaL + Ici{Ca) + IsDc) 



Ixir and determine the membrane potential of the cell at rest near -70 mV and are 
specified in 

The equation describing the L-type Ca-current {Icai) in terms of the Hodgkin-Huxley 
kinetics of the L-type Ca-channel, is given by 

IcaL =mhVCa GcaL{Vm " Ecal), (2) 

where m is the voltage-dependent activation variable, h is the volt age- dependent inactivation 
variable and Vca is the inactivation parameter. The dynamics of the variables m and h are 
described by first order differential equations of the Hodgkin-Huxley type The calcium- 
dependent inactivation is given by vca = K^cal i\Ca^cyt\ + K^ca)- 
The Ca-dependent Cl-current Ici(Ca) is given by 



Ici{Ca) - r„ 2+1 , "^t} Gci(Ca) (Ki - Eci(Ca)) (3) 

The chloride current increases with the cytosolic calcium concentration [Ca^^j^], causing a de- 
polarization to the Nernst potential of chloride ions {Eci{Ca)) near -20 mV in NRK fibroblasts 
for sufficiently high values of [Ca^jj^]. 

The store-dependent calcium current Isdc is described by 

hoc = r„ 2+/!^^^ GsDc{Vm- EsDc)- (4) 

[Ca^+ ] + KsDC 

This store-dependent calcium channel allows calcium ions to fiow from the extracellular s pac e 
into the cytosol at a rate inversely proportional to the calcium concentration in the ER ()28| ). 
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SDC channels are thought to play a major role in the control of Ca-homeostasis in the cell 
(22, 3). 



The rate of change of Ca-content of the cytosol of the cell due to inflow through the cell 
membrane and from the ER store, and by buffering is described by 

yolcyt——^^^ = ApM JpM + Aer{,JiPzR + JlkER — JsERCa) — Volcyt — ^' (^) 

where Volcyt represents the cytoplasmic volume and ApM and Aer the area of the cell 
membrane and of the ER membrane, respectively. The term [BCa] denotes the buffer- 
calcium complex in the cytosol and will be explained later. The flux of calcium through the 
membrane [Jpm) is the sum of the influxes of Ca^"*" ions through the L-type Ca-channel, 
through the SDC-channel, and of the extrusion by the PMCA-pump j^), and is given by 

JpM = —{i/{zcaFApM)){IcaL + IsDc) — J PMC A- 

The djTiamics for the intracellular calcium oscillator is described by the flux of calcium 
through the ER membrane. The rate of change of calcium content in the ER depends on 
the sum of flux through the IPs-receptor {Jip^p), flux by leak through the ER-membrane 
{JlkER) and flux by removal by the SERCA pump {Jserca), which results in 

VoIer ^ = Aer{—JipzR — JlkER + Jserca), (6) 

where VoIer represents the volume of the ER. 

The flux through the IPs-receptor is described by 

JmR = foo Kjp,n {[Cal\] - [Cali]) (7) 

where [Ca^] — [Ca^jjt] is the concentration difference between calcium in the ER and in the 
cytosol. Kjp^R is the rate constant per unit area of IPs-receptor mediated release. The terms 
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/oo and w represent the fraction of open activation and inactivation gates, respectively, /o, 
and Woo depend both on the cytosohc calcium concentration and are described by 



KfiP, + [Cai+] 
and 



""^ [IPs] , ^ \rn^+] ■ ^ ' 



The inactivation time constant of the IPa-receptor is defined by 



(10) 



Kfjp^, Kyjip^, Ky^i^Ca) and a are constants. The fraction of open activation gates (f) is inde- 
pendent of the IP3 concentration, but increases when the calcium concentration in the cytosol 
increases. The fraction of open inactivation gates (w) depends on the IP3 concentration and 
on [Ca^jjf]. Tw determines the duration of the de- inactivation of w. 

JikER is a passive leak of Ca^+ from the ER into the cytosol which is not mediated 
by the IPs-receptor, but by an additional Ca-channel in the ER membrane, presumably the 
translocon. Experimental evidence for a role of the translocon complex as a passive Co?^ leak 
channel has been presented recently JikER is given by JikER = KikER{[Co?j^p\ — [Ca^y^]). 
We used the leakage parameter KikER as a control parameter to study the dynamics of the 
plasma membrane, because changes in the leak of Ca-ions through the ER membrane produce 
proportional changes in [Ca^^]. 

JsERCA represents the flux of calcium into the ER by the SERCA pump and is given by 

JsERCA = JsERCA {[Calytf / {K'^erca + [Calytf)}. 

Finally, calcium in the cytosol is buffered by proteins in the cytosol. The dynamics of 
buffering is given by d[BCa]/dt = kon{[TB] - [BCa])[Cal+t] - Kff [BCa], where [Tb] is the 



8 



total concentration of buffer in the cytosol and fcon and k^fj are the buffer rates 

The excitable membrane and the IPs-mediated intracellular calcium oscillator are cou- 
pled by the Ca-concentration [Calyf] in the cytosol as explained in (!22i ). During an action 
potential, opening of the L-type Ca-channel causes a large inward current of Ca-ions through 
the plasma membrane. The increased [Ca^^j^] activates the IPs-receptor (calcium release 
channel), causing calcium release from the ER, which further contributes to the intracel- 
lular cytosolic calcium transient. In the reverse process, IPs-mediated calcium oscillations 
cause periodic calcium transients, which lead to periodic opening of the Ca-dependent Cl- 
channels. The depolarization of the membrane potential towards the Nernst potential of 
the Ca-dependent Cl-channels near -20 mV causes activation of the L-type Ca-channels in 
the plasma membrane and excitation I31I ). After an action potential or Ca-transient 
the reduction of cytosolic calcium by the activity of the SERCA and PMCA pumps reduces 
Ici{Ca) (see Eq. IH]), enough to allow the membrane to return to the membrane potential at 
rest near —70 mV. 

The dynamics of the single-cell model depends on seven variables (m, h, w, [BCa], V^, 
[Caly-i-] and [Ca^"^]), which were defined above. To study the stability of the complete system 
we have determined the sin gula r states for the system and calculated the Floquet multipliers 



of these singular states (j32 



331). 



Stability analysis of the membrane model 

We will first analyze the bifurcations and local stability of both the excitable membrane and 
intracellular calcium oscillator separately, and then compare the results with the properties 
of the single-cell model including both the membrane dynamics and intracellular calcium 
oscillator. Two different analyses, namely, our own implementation in C, and the software 
package XPPAUT (34), which includes an AUT086 (2^ interface, gave the same results. 

In the single-cell model the intracellular calcium oscillations can be eliminated by setting 
the IPs concentration ([IPs]) to zero. This allows the study of the excitable cell membrane 
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separately from the calcium oscillator. The dynamics of the plasma membrane depends 
on the cytosolic calcium concentration [Ca^^j^], which opens the Ca-dependent Cl-channel. 
Since the leak of Ca-ions from the ER affects the mean value of [Ca^jj^], the dynamics of the 
membrane is studied as a function of the leakage parameter Kiuer- Fig. |21 shows a hysteresis 
diagram for the excitable cell membrane with the steady states of the calcium concentration 
in the cytosol ([Ca^^j^], panel A) and of the membrane potential (Kn, panel B). The thick 
and thin solid lines refer to the stable states for increasing and decreasing values of KikER, 
respectively. The dashed-dotted hues reflect the transitions between the two stable branches 
for increasing and decreasing values of Ki^er- 

Starting at the value zero for Ki^er, the inwardly rectifying K-channels keep the mem- 
brane potential at the resting membrane potential of the NRK fibroblasts near —70 mV, 
where the membrane is able to produce an action potential upon electrical stimulation 
For increasing values of Ki^er, [Co^jj^] and Vm increase gradually, causing a decreasing thresh- 
old for activation. The gradual increase of Vm is due to gradual opening of the Ca-dependent 
Cl-channels for increasing [Ca^^] (see Eq. E)). At Ki^er ~ 58.0 x 10^^ dm/ s, [Ca1yf\ is large 
enough to open the Ca-dependent Cl-channels driving the membrane potential towards the 
Nernst potential for Cl~-ions which is near —20 mV in NRK fibroblasts (see Fig. I2l3). The 
resulting depolarization causes closure of the inwardly rectifying K-channels and opening 
of the L-type Ca-channels which leads to an increase of calcium inflow from the extracel- 
lular medium into the cytosol. The positive feedback via the membrane potential between 
Ca-dependent Cl-channels and L-type Ca-channels explains the abrupt increase of [Ca^jj^] 
(dashed-dotted line) to 2.3 /iM. 

When we decrease Ki^er starting from 60.0 x 10^^ dm/ s (thin solid line), the cell remains 
depolarized near -20 mV far below the value of KikER at 58.0 x 10~* dm/ s. This is caused 
by the feedback between the Ca-dependent Cl-channels and L-type Ca-channels. When 
KikER decreases, [Ca^jji] also decreases, which reduces the fraction of open Ca-dependent 
Cl-channels. As a consequence, the membrane potential slightly decreases just below -20 
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mV, which leads to an increased fraction of open L-type Ca-channels, since the product of 
mh of steady-state activation and inactivation (see Eq. 12)) reaches a maximum just below 
-20 mV. The increment of the fraction of open L-type Ca-channels leads to an extra inflow 
of calcium in the cytosol, which increases the fraction of the open Ca-dependent Cl-channels 
and prevents the system from falling back to a membrane potential near —70 mV. Thus, 
in spite of the slow decrease of calcium concentration and membrane potential caused by 
the leak channels, the feedback by the L-type Ca-channels keeps the system at an elevated 
[Ca^^] and membrane potential near —20 mV until low values KikER- The calcium in the 
cytosol returns to a low concentration, only when the Ki^er is decreased to very low values. 
Then the Ca-dependent Cl-channels close and the membrane potential repolarizes to —70 
mV. Separate simulations showed that the inward rectifier contributes to the transitions, 
but not to the hysteresis. 

Stability analysis of the intracellular calcium oscillator 

Following a similar plan as for the excitable cell membrane, we obtained a bifurcation diagram 
for the intracellular calcium oscillator as a function of the IP3 concentration under conditions 
that the L-type Ca-channels are blocked such as with nifedipine. This was achieved by 
setting GcaL to zero and KikER to its physiological value of 2.0 x 10^^ dm/s. In this way, 
we eliminate the contribution of calcium inflow by the L-type Ca-channels and remove a 
principal influence of the membrane model on the intracellular calcium oscillator. Therefore, 
we only take into account the Ca-flux through the SDC-channels and PMCA pump in the 
plasma membrane. 

AS e.p.a.ed . Q, .>..e of PMCA ... SEKCA p„.„p . c™> 

to reproduce the steady state calcium concentrations in the cytosol and in the ER. By 
eliminating the calcium inflow by the L-type Ca-channels, less calcium flows into the cell. 
Therefore, we have to change the relative strength of the PMCA and/or SERCA pump to 
maintain the proper balance between calcium concentration in the cytosol and ER. In this 
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model study, we choose to decrease the strength of the SERCA pump. By doing so, the 
system reveals a bifurcation diagram (Fig. EK) similar to that observed in other models 

0,0). 

Fig. shows the dynamical behavior of [Ca^^j^] as a function of [IP3] for Gcai = 
and with Js^erca set to 2 x 10^^ {fimol)/{s x drn^). Fig. ElA. shows a single stable steady 
state for small values of [IP3] (range 0.0 - 0.2 /iM). At [IP3] near 0.2 /xM the dynamics 
reveals a supercritical Hopf bifurcation (thick solid line), and the system becomes a calcium 
oscillator in the range of IP3 concentrations between 0.2 and 3.6 /xM. For [IP3] values near 
3.6 /iM the system meets a supercritical Hopf bifurcation and remains stable for higher 
IP3 concentration at a Ca-concentration near 5 /iM. In the range for [IP3] above 3.5 /iM, 
the elevated mean level of [Ca^^j^] gives rise to a short time constant for the inactivation 
parameter w (Eq. [TUI) . Due to this small time constant the inactivation w recovers relatively 
fast compared to the removal of [Ca^^], i.e. before the activation parameter / de-activates 
to small values. As a result the product fw does not reach small values and the IP3-receptor 
remains open (see Eq. [7j), causing a constant leak of calcium. 

For decreasing [IP3] values (thin solid line), the system starts at a stable fixed point 
which remains stable until 3.50 /iM. In the range [IP3] between 3.5 and 3.6 /iM, the system 
exhibits bistability and a hysteresis over a small range of IPs-values. This hysteresis is caused 
by the positive feedback between [Ca^jj^] and the activation gate (/). For decreasing [IP3], 
the [Ca^^] is already elevated and so a large fraction of activation gates / is already open and 
the time constant r^j is short. Due to the short time constant r^, the time for de-inactivation 
(w) is faster than for de-activation (f). As a consequence, the product of / and w does not 
reach small values and calcium passes continuously through the IPs-receptor from the store 
into the cytosol. This hysteresis did not show up in the figures presented by Li & Rinzel 



f)36[ ). but appears in their model if we insert the parameter values which apply to the NRK 
fibroblasts (see ([3)). 

For [IP3] values below 3.5 /iM, the stable fixed point disappears, and the system starts 
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to operate as an oscillator, until [IP3] values smaller than 0.15 /iM, where the system returns 
to a single stable steady state. 

As a next step, we have set the strength of the SERCA pump back to its default value 
8.10"^ fimol/{s X dm'^) which corresponds to the value in the single-cell model with an 
excitable membrane and IPs-mediated calcium oscillations. This results in the bifurcation 
diagram shown in Fig. |2|3. Fig. |2P shows a major hysteresis in the [IP3] range between 8 
and 53 /iM (see inset). To compare the results with those in Fig. we scaled Fig. in 
the same [IP3] range as in Fig. IHK. 

When the strength of the SERCA pump is increased to 8.10^^ yumo//(s x dm^), [Caly^] 
decreases more rapidly after a calcium transient. This affects the time constant of the 
inactivation parameter w (see Eq. HHI). For small [Ca^+J levels, this time constant is relatively 
large, ensuring a slow de-inactivation. This explains why a more powerful SERCA pump 
gives rise to calcium oscillations over a much larger range of IP3 concentrations. Only at 
sufficiently large [IP3] values does become sufficiently small such that de-inactivation (w) 
takes place more rapidly than de-activation (/). For these high IP3- values, the product fw 
of the activation parameter (f) and the inactivation parameter (w) is large enough to allow 
a continuous leak of calcium through the IPs-receptor. 

The inset in Fig. 133 shows a single stable steady state for small values of [IP3]. At [IPs] 
near 0.2 fiM, the dynamics reveals a subcritical Hopf bifurcation (thick solid line), and the 
system becomes a calcium oscillator in the range of IPs concentrations between 0.2 and 53 
fiM. For [IPs] above 53 fiM, [Ca^+] is elevated at a steady state concentration near 4 /xM. 
For decreasing [IP3] values (thin solid line), the system starts at a stable elevated [Ca^^] 
which remains stable until [IPs] is near 8 fiM. For [IPs] values below 8 fiM, the stable fixed 
point disappears, and the system starts to operate as an oscillator, until [IPs] values smaller 
than 0.2 fiM. We conclude that increasing the activity of the SERCA pump makes it more 
easy for the cell to oscillate at higher [IPs] values and causes a hysteresis over a larger range 
of [IPs] values. 
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Stability analysis of the single-cell model 

Unblocking the L-type Ca-channels {GcaL =0.7 ^M) transforms the bifurcation diagram of 
EfB into that of Fig. Fig. |3] shows [Ca^^j^] (panel A) and the membrane potential (panel 
B) as a function of IP3 concentration in the cell. The solid and dashed-dotted lines represent 
stable and unstable states, respectively. For small [IP3] values in the range from 0.00 to 0.15 
/iM, the cell has a single stable steady state ("resting state") with a membrane potential near 
-70 mV. For [/P3] > 0.15 /iM, the stable fixed point becomes unstable in a subcritical Hopf 
bifurcation. Calcium oscillations together with action potentials occur for IP3 concentrations 
in the range between 0.15 and 1.75 /iM ("AP-firing state") (see panel C which shows the 
membrane potential as a function of time for [IP3] = 0.7 /iM). In this regime, a rapid calcium 
inflow from the ER into the cytosol opens the Ca-dependent Cl-channel, causing an inward 
current towards the Cl-Nernst potential close to —20 mV . This depolarization activates the 
L-type Ca-channels leading to an AP. After closure of the IP3-receptor, calcium is removed 
from the cytosol by the Ca-pumps in the cell membrane and ER, leading to repolarization to 
-70 mV. For [IP3] > 1.75 /xM, the fixed point ([Ca^+], V^) near (3.00 /iM,-20 mV) becomes 
stable in a subcritical Hopf bifurcation ("depolarized state"). This can be understood from 
the fact that the time-constant (see Eq. ITn|l for calcium-dependent (de-)inactivation 
of the IP3-receptor decreases for increasing values of [IP3] and for increasing values of the 
mean [Ca^^j^]. Near [IP3] = 1.75, the time-constant is relatively short. During a cytosolic 
calcium transient the fast inactivation of the inactivation gates w of the IP3-receptor is 
followed by a fast de-inactivation of the inactivation gates of the IP3-receptor. During the 
fast de-inactivation, the fraction of open activation gates / of the IP3-receptor is still high 
due to high [Ca^^j^] (because removal of calcium through the SERCA and PMCA pump is 
not fast enough). As a consequence the IP3-receptor remains open. Now the /P3-receptor 
acts as a constant leak channel, like Jiuer- This leak of calcium into the cytosol opens the 
Ca-dependent Cl-channels, causing a maintained depolarization to the Cl-Nernst potential 
near —20 mV ("depolarized state") (panel B). 
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If [IP3] is decreased starting from [IP3] =2.5 /iM, the cell with both the excitable 
membrane and intracellular calcium oscillator active exhibits a complex hysteresis pattern. 
For decreasing IP3 concentrations, the system stays in a single stable state (solid line) at 
an elevated [Ca^^] near 3 /iM and a membrane potential near —20 mV until [IP3] ^ 0.85 
fiM ("depolarized state"). Then, the cell goes through a Hopf bifurcation (dashed line) 
forcing the system to behave as a stable oscillator with small calcium oscillations with an 
amplitude of about 6 /iM and with small membrane potential oscillations around —23 mV 
("fast oscillating state"). These small oscillations of the membrane potential just below —20 
mV are illustrated in more detail in panel D. Note that the oscillations of the membrane 
potential in panels C and D are both obtained for [IP3] =0.7 /iM, illustrating the hysteresis. 
The oscillations shown in panel D are due to small IPs-mediated calcium oscillations with 
active involvement of the L-type Ca-channel dynamics. Due to decreasing [IP3], the de- 
inactivation time constant of the IPs-receptor increases gradually. This makes it possible 
for the cell to generate calcium oscillations. Setting Vca in the equation for the L-type Ca- 
current to 1 does not change the bifurcation diagram of Fig. |3] The shape of the bifurcation 
scheme in Fig. 1314. remains the same, but the calcium oscillations extend over a larger range 
of [Ca^jjj]-values (approximately twice as large). 

At [IP3] ~ 0.45 yuM the stable small-amplitude oscillator becomes unstable (dashed line), 
returning the system to the stable oscillations with large amplitude Ca-oscillations with a 
peak value near 20 fiM and with action potentials in the range between -70 and -10 mV 
("AP-firing state"). Finally, for [IP3] values smaller than 0.15 fiM the system returns to a 
single stable state ("resting state"). 

In comparison with the simple dynamics of the cell membrane and intracellular Ca- 
oscillator, shown in Figs. El and El it is remarkable to see the complex behavior of the 
single-cell model shown in Fig. |31 

Since the SDC channels in the plasma membrane play a crucial role in stabilization of 
...... fl Q, we s..e. . ceU . a ...«o„ o. SDC 
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conductance in a range between 0.00 and 0.20 nS. Fig. |31 shows the hysteresis diagrams for 
five different values of Gsdc- As explained in the calcium homeostasis of the cell is 
unstable for Gsdc = 0.00 nS. For small values of Gsdc bistability and hysteresis appears. 
The IP3 range with hysteresis is largest for a Gsdc value near 0.04 nS (see Fig. El). For 
higher values of Gsdc, the range of hysteresis decreases until the typical Hopf-bifurcation 
for the intracellular IPs-mediated calcium oscillations remains for Gsdc = 0.20 nS. 

The IP3 range of the hysteresis as a function on the SDC conductance channel is shown 
in Fig. ini We define the IP3 range of hysteresis as the range of [IP3] in which multiple states 
are found for increasing and decreasing [IP3]. For example, in Fig. |3 hysteresis takes place 
for [IP3] values between 0.45 and 1.75 fiM, giving an IP3 range of hysteresis of 1.3 fiM. 

Recent data in the literature show that the SDC conductance, which was found to give 
the largest range for hysteresis in our study (near 0.04 nS), corresponds to the observed SDC 

41). The SDC conductance reported in Table 1 of (js^ 



conductance in other studies 
and in 



HQ 

( 401 ) was in the range between 0.04 and 0.05 nS (solid line below the peak in Fig. 



in)). For hepatocytes ()4lh a SDC conductance was reported in the range between 0.08 and 



0.14 nS. However, since the density of all ion channels in hepatocytes is twice as high as in 
fibroblast (42, 3), the ratio of conductances for the ion channels is the same in hepatocytes 
and NRK fibroblasts. If we correct for this higher density, rescaling all conductances for 
those of NRK fibroblasts, we obtain the dotted line in IHl Therefore, the SDC conductance, 
for which hysteresis is found over the largest range of IP3- values in our study (see Figs. El 
andEI), is in agreement with experimental observations for SDC conductance. 

Discussion 

In this study we have analyzed a relatively simple model with an excitable membrane and 
with IPs-mediated calcium oscillations. The interaction between these mechanisms in a 
single-cell model revealed a surprisingly rich behavior with stable/instable states with hys- 
teresis. 
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The hysteresis and bistabihty of the membrane potential and the intracellular calcium 
concentration in Fig. ^A. and B obtained by stability analysis of the single-cell model provide 
an explanation for the various growth-state dependent changes in the electrophysiological 
behavior of normal rat kidney (NRK) fibroblasts in cell culture (3)- The stability analysis 
of the single-cell model (Fig. reveals for low [IP3] values (range 0.0 - 0.2 fiM) a cell in 
the "resting state" . Increasing the [IP3] leads to spontaneous AP firing ("AP-firing state") 
and at high [IP3] values above 1.75 fiM, the cell depolarizes ("depolarized state"). When 
we start at an [IP3] value of 2.0 /iM and decrease [IP3], the system is in the "depolarized 
state" and changes from the "fast oscillating state" (range 0.45 - 0.8 fiM) to the "AP-firing 
state" (range 0.2 - 0.45 fiM) back to its "resting state" ( rang e 0.0 - 0.2 fiM), which is in 



agreement with experimental data shown by Harks et al. ()27l ). In the study of Harks et al. 
( 23^ washout of the medium conditioned by the transformed cells by perfusion with fresh 
serum-free medium, causes the cells to slowly repolarize, and, preceded by a short period of 
fast small-amplitude spiking of their membrane potential ("fast oscillating state"), to regain 
spontaneous repetitive action potential firing activity similar to that of the contact inhibited 
cells. This compares well with the results in Fig. 1^3, which shows for decreasing [IP3] a very 
similar behavior. Therefore, we conclude that the stable states of the model as revealed by 
stability analysis of the single-cell model correspond in great detail to the observed growth- 
state dependent modulations of the membrane potential of NRK fibroblasts in cell culture. 
This strongly suggests that these growth-state dependent modulations of the membrane 
potential of NRK cells reflect just different states of the same cell, rather than the behavior 
of cells, that have differentiated to different cell types with different properties during the 
various stages of growth-factor stimulated development in vitro. 

I I 

Most of the parameter values in our model were taken from the literature (see ff22>) ) for 
a detailed overview). Interestingly, the parameter values for the excitable membrane and for 
the IP3-mediated calcium oscillator, which are very different mechanisms, are not indepen- 
dent. This can be understood from the fact that the dynamics of the excitable membrane 
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and of the IP3 receptor are coupled by the cytosohc calcium concentration. Changing one 
parameter of the excitable membrane or calcium oscillator affects the other mechanism by 
changes in the cytosolic calcium concentration. This is illustrated, for example, by Fig. IHl 
Changing the strength of the SERCA pump causes large differences in the range of hys- 
teresis in the dynamics of cytosolic calcium (Fig. E)), and therefore also in the dynamics 
of the membrane potential (see Fig. HJ which illustrates the relation between the dynamics 
of cytosolic calcium and the membrane potential). Although the parameter values for the 
excitable membrane and for the IPs-mediated calcium oscillator were taken from different 
studies, they fit nicely together to explain the behavior of NRK cells both qualitatively and 
quantitatively. This provides strong evidence for the reliability of these parameter values. 
Moreover, this suggests that cells should have complicated regulatory mechanisms to control 
all parameter values within a proper range of parameter values to ensure the proper cell 
dynamics. 

Summarizing, we explored the dynamical properties of a single-cell model reproducing 
experimental observations on calcium oscillations and action potential generation in NRK 
fibroblasts. A bifurcation analysis revealed hysteresis and a complex spectrum of stable and 
unstable states, which allows the system to switch among different stable branches. Sta- 
bility of the cell behavior is dominated by the homeostatic function of the SDC channel. 



The conductance, which provides the largest 
experimental values for this conductance 



Pa range for hysteresis, compares well with 



m 



4ll ). Experimental observations in NRK 



fibroblasts revealed the same kind of hysteresis as shown by this study. 
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Figure Legends 



Fig.m 

Conceptual model of the processes involved in membrane excitability and intracellular IPs- 
mediated calcium oscillations by calcium release through IPs-receptors in the ER mem- 
brane in NRK fibroblasts. Cell-membrane excitability is supported by inwardly rectifying 
potassium channels (GKir), Ca-dependent Cl-channels {Gci{ca)), L-type Ca-channels (Gcai), 
store-dependent calcium (SDC) channels (Gsdc), a PMCA pump and leak channels (Gik)- 
The total fiux of calcium through the ER-membrane is the result of the contribution by 
the SERCA pump, by the IPs-receptor (J/pgi?) and by leak channels in the ER membrane 
{JikERj- The membrane excitability and IPs-mediated calcium oscillations are coupled by 
the cytosolic calcium concentration, which is also affected by a calcium buffer B. 

Fig. El 

Stable and unstable states for the excitable membrane using KikER as a control parameter. 
The intracellular Ca-oscillator was silenced by setting the IPs concentration to zero. Thick 
(thin) lines correspond to the stable steady-state solutions for the [Ca^^] (panel A) and for 
.he .e.b.a„e poten.M (p^e, B, fo. ,n.ea.n. (dec....) values of fo. .... The se. 
of parameter values in this model was as reported in (j22i), with Gsdc = Gsoc ■ 

Fig. [31 

The bifurcation diagram for the intracellular calcium oscillator in the single-cell model as a 
function of [IPs] after elimination of action potentials {GcaL = 0) and with Js^erca set to 2 
X 10^^ fimol/{s X drn^) in panel A and to 8 x 10^^ fimol/{s x drn^) in panel B. Analogous to 
Fig. 121 thick and thin solid lines correspond to the stable states for increasing and decreasing 
values for [IPs], respectively. The three insets in panel A show the Ca-concentration as a 
function of time for [IPs] values at 0.01, 2 and 4 /iM. The inset in B shows the stable (solid 
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lines) and unstable (dashed-dotted lines) states for a large r ang e of [IP3] values. The set of 



all other parameter values in this model was as reported in 



22), with GsDc — Gsoc- 



Fig. m 

The bifurcation diagram for the single-cell model. The figure shows the stable (solid lines) 
and unstable (dashed-dotted lines) states for [Ca^jj^] (panel A) and the membrane potential 
(panel B) as a function of IP3 concentration. Panel C and D show the membrane potential 
for [IP3] at 0.7 /iM in case of increasing and decreasing [IP3], respectively. The small arrows 
on the curves show the direction of change of the stable modes for increasing and decreasing 
values of [/-P3]. The set of parameter values in this model was as reported in js^. 



Fig.il 

Bifurcation diagrams for the single-cell model as shown in Fig. 01^. for different values of the 
SDC channel conductance from Gsdc = 0.02 nS (bottom) to Gsdc = 0.20 nS (top). At 
0.04 nS (second graph from bottom), the hysteresis loop has a maximum in the [IP3] range 
from 0.5 to 1.95 /iM. 

Fig. El 

IP3 range of the hysteresis loop in /xM (cf. Fig. E)) as a function on SDC conductance (Gsdc) 
for a NRK cell with a capacitance of 20 pF reveals a value in which the hysteresis area is 
maximum (dots are simulated data). Below the peak, we illustrate the SDG conductance 

nil n 

measured in experiments and separately reported in ()39l 1401 ) (solid line) and ()41f ) (dashed 
line) . 
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Figure 4: 
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